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The expected amplitudes and spectral slopes of relic gravitational waves, plus their specific 
I i correlation properties associated with the phenomenon of squeezing, may allow the registra- 

^ • tion of relic (squeezed) gravitational waves by the first generation of sensitive gravity-wave 

Q\ ' detectors. 
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Relic gravitational waves should have been produced by strong variable gravitational field 
of the very early Universe which parametrically (superadiabatically) amplified the zero-point 
. quantum oscillations of the gravitational waves [1] . The generating mechanism is universal and 

*q | relies only on the validity of the general relativity and basic principles of quantum field theory. 

The initial vacuum quantum state of each pair of waves with oppositely directed momenta has 
evolved into a highly correlated multiparticle state known as the two-mode squeezed vacuum 
quantum state [2] . (For a recent review of squeezed states see [3] ) . The phenomenon of squeezing 
manifests itself now in a specific standing-wave pattern and periodic correlation functions of the 
^ generated field [4] . It is remarkable that this statistical signature may significantly facilitate the 

detection of the relic (squeezed) gravitational wave background. It is possible, as we will show 
below, that the appropriate data processing will allow detection of relic gravitational waves by 
first generation of the forthcoming sensitive instruments, such as initial laser interferometers in 
LIGO [5], VIRGO [6], GEO600 [7]. 

We consider cosmological gravitational wave field h%j defined by the expression 

ds 2 = a 2 {r])[dr, 2 - (<% + h ,j)d.r\\.r J . (1) 

The Heisenberg operator for the quantized real field hij can be written as 



L(??)e inx c n + C(r?)e- inx c n 



(2) 



where C = \A67r Ipi and the Planck length is lp\ = (Gh/c 3 ) 1 ^ 2 , the creation and annihilation 

operators satisfy c n |0) = 0, [c n , c m ] = 5 ss '5 3 (n — m), and the wave number n is related to the 
wave vector n by n = (Sijtfni) 1 / 2 . The numerical value of the constant C is determined by the 
quantum normalization, that is, by the requirement that, initially, each mode of the field has 

s 

had the energy of only of a "half of the quantum". The Py(n) (s = 1,2) are two polarisation 

s 

tensors. The functions h n {v) are governed by the linearised version of the Einstein equations. 



For every wave number n and each polarisation component s, the functions hn(v) have the form 

® s s ^ s s 

/in(^) = [ n n(f?) + v n{v)]/ a (v)i where u n (r)) and v n (r]) are expressed in terms of the three real 
functions (for each s): r n - squeeze parameter, 4> n - squeeze angle, n - rotation angle, 

u n = e Wn cosh r n , v n = e ~ i( - 9n ~ 24 '^ sinh r n . 

The functions r n (rj), 4> n (ri), d n (rj) obey the dynamical equations [4]: 

r' n = — cos 2(f) n , (f>' n = —n — — sin 2<p n coth 2r n , 9' n = — n — — sin 2(p n tanh r n , (3) 

where ' = d/drj, and the evolution begins from r n = 0, which characterizes the initial vac- 
uum state. These equations can be transformed into a more familiar form of the second-order 

S s s* 

differential equation for the function M n (??) = u n (rj) + v n {r]) [1]: 

= 0. (4) 

Dynamical equations and their solutions are identical for both polarisation components. 

The present day values of r n and (p n are essentially all we need to calculate. The mean number 
of particles in a two-mode squeezed state is 2 sinh 2 r n (for each s). This number determines the 
mean square amplitude of the gravitational wave field. The time behaviour of the squeeze angle 
4> n determines the time dependence of the correlation functions of the field. The amplification 
(that is, the growth of r n ) governed by (3) is different for different wave numbers n. Therefore, 
the present day results depend on the present day frequency v {y = cn/2ita) measured in Hz. 

In the long-wavelength regime, that is, during the interval of time when the wavelength 
A (77) = 2ira/n is longer than the Hubble radius l(rf) = a 2 /a', the squeeze parameter r n (rf) 
grows with time according to r n (rj) pa log [a(r])/a*], where a* is the value of a(rj) when the long- 
wavelength regime, for a given n, begins. The final amount of r n is r n pa log [a** /a*], where a** is 
the value of a{jf) when the long-wavelength regime and amplification come to the end. (During 

s 

the long- wavelength regime, the dominant solution for the function hn(v) i s approximately 
constant, instead of the adiabatic decay h n {v) ~ e a (2\ which takes place in the short-wavelength 
regime.) After the end of amplification, the accumulated (and typically large) squeeze parameter 
r n stays approximately constant, and the functions ^> n (ry), 9 n (rj) are 4> n = —n(r] + rj n ), 8 n = 4> n . 
(For more detail see [8] and references therein.) The constant % is determined by the moment 
of time when the long-wavelength regime, for a given n, terminates. This constant varies from 
one wave number n to another, but does not change too much at the intervals An pa n. This 
constant will eventually survive as the (not very important) additive constant t v in phases of 
the presently existing oscillating field cos[27rz/(i — t v )], where t is the total elapsed time (billions 
of years since the end of the amplifying regime for the mHz — kHz frequencies of the current 
experimental interest), and therefore v(t — t v ) >> 1 for these frequencies. 

The numerical results depend on the concrete temporal behaviour of the gravitational pump 
field represented by the cosmological scale factor a{rj). We know that the present matter- 
dominated stage a(rj) oc rj 1 was preceeded by the radiation-dominated stage a{rj) oc rj. Both 
these functions are power-law functions in terms of rj time. The function a{rj) describing the 
initial stage of expansion of the very early Universe (before the era of primordial nucleosynthesis) 
is not known. It is convenient to parameterize the a(rj) at the initial stage also by power-law 
functions of 77. It is known [1] that power-law functions a(rj) produce gravitational waves with 
power-law spectra in terms of u, and, vice versa, every piece of the generated spectrum which 
can be approximated by a power-law function in terms of v has been generated by a piece of 
evolution which can be approximated by a corresponding power-law function a(rj). 
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Concretely, we take 0^(77) at the initial stage of expansion as 0,(77) = ^oM 1+/3 where 77 time 
grows from —00. Expansion, that is increase of a(rj), at this interval of evolution requires the 
parameter (3 to be (3 < — 1. This parameter will eventually determine the spectral slope of 
the generated spectrum. The parameter l Q has the dimensionality of length and is associated 
with the Hubble radius (in general, time-dependent) at the initial stage of expansion. In the 
special case (3 = —2 one encounters a portion of the De Sitter (inflation) evolution, whereby 
the Hubble radius remains strictly constant and equal to l . The ratio lpi/l characterizes the 
"strength" of the pump field and will eventually determine the mean square amplitudes of the 
generated field. From 77 = 771 , 771 < 0, the initial stage of expansion is followed by the radiation- 
dominated stage a e {rj) = l a e (ri — i] e ) and then, from 77 = 772 , by the matter-dominated stage 
— Tim) . The constants a e , a m , rj e , r/ m are expressed in terms of the fundamental 
parameters l Q , (3 through the continuous joining of a{rj) and a'(rf) at 771, 772- The present era is 
defined by the observationally known value of the Hubble radius ljj = c/i/~ 2 x 10 28 cm. We 
denote this time by T]n and choose t]r — i] m = 1, so that <i(t]r) = 21h- The ratio a(r]R) / '0(^2) = z 
is believed to be around z = 10 4 . The wave numbers n#(n# = 47r), n m (n m = y/znu), n c denote 
the waves which are leaving the long- wavelength regime at, correspondingly, rjR, 772, 7/1 • The 
shorter waves, with n > n c , have never been in the amplifying long-wavelength regime. (The 
present day frequency is c , corresponding to the wave number n c , is around 10 10 Hz.) 

The most complete description, allowed by quantum mechanics, of the generated field is 
provided by the evolved operator ^(77, x), eq. (2). The field is stochastic in the sense that it is 
characterised by the quantum mechanical mean values, variances, various correlation functions, 
etc.. The mean value of the field hij is zero, (0\hij(rj, x)|0) = 0, at every moment of time and in 
each spatial point. The variance (0|/ijj(?7, x)^*- 7 (77, x)|0) = (h 2 ) is not zero, and it determines 
the mean square amplitude of the generated field - the quantity of interest for the experiment. 
Taking the product of two expressions (2), one can show that 

(h 2 )= rV(n>i7)-, 
Jo n 

where, for the present era, 

f n \ l3+2 

h(n) « A I , n < tih, (5) 

\nH/ 

/ n \P 

h(n, 77) w A cos 4> n (77) — , n H <n<n m , (6) 
\riHj 

( n fnu\ 

h(n, 77) « A cos <fi n (v) [ — — , n m < n < n c , (7) 

\n H J \n m J 

and 

Ipi 8v^+ 2 



A = 



io li + z?!^ 1 ' 

The periodic time-dependent structure cos n (ry) is a consequence of squeezing, and it will be 
analyzed below. 

So far, the parameters l and (3 were quite arbitrary. The available information on the 
microwave background anisotropies [9, 10] allows us to obtain some information about A and (3. 
The quadrupole anisotropy produced by the spectrum (5) - (7) is mainly accounted for by the 
wave numbers near njj. Thus, the numerical value of the quadrupole anisotropy produced by 
relic gravitational waves is approximately equal to A. Since (according to [11]) the quadrupole 
contribution of relic gravitational waves is not smaller than that produced by primordial density 
perturbations, this gives us A pa 1CP 5 . We are not certain whether a significant part of the 
quadrupole signal is indeed provided by relic gravitational waves, but we can assume this. The 



evaluation of the spectral index n of the primordial perturbations resulted in n = 1.2 ± 0.3 [10] 
or even in a significantly higher value n = 1.84±0.29 [12] (see also [13], where one of the best fits 
corresponds to n = 1.4). The most recent analysis [14] of the observational data favors n = 1.2 
and the quadrupole contribution of gravitational waves twice as large as the contribution of 
density perturbations. We interpret these evaluations as an indication that the true value of n 
lies somewhere in the interval n = 1.2 ~ 1.4 (hopefully, the planned new observational missions 
will determine this index more accurately). Since the primordial (before further processing in 
short-wavelength regimes) spectral index n is related with (3 by n = 2(3 + 5, this gives us the 
parameter (3 in the interval (3 = —1.9 ~ —1.8. The derived parameters A, (3 define the numerical 
level of h(n,rj) in the frequency domain accessible to laser interferometers (eq. (7)) and make 
the gravitatinal wave signal measurable, as we will discuss shortly. 

One can often hear that "inflation predicts" a negligibly small contribution of gravitational 
waves, and, as a result, hopeless prospects for direct detection of relic gravitational waves. This 
is a damaging statement, and it requires special attention. The origin of this "conventional 
wisdom" is the (incorrect) prediction of inflationary theorists for the amplitudes of density 
perturbations, which, in addition to gravitational waves, can also be produced, under some 
extra assumptions, by the same amplifying mechanism. The central inflationary formula for the 
generated matter density variation (or curvature perturbation) is 

Sp 1/3/2 

7 ~ T ~ ~v r ' ( ) 

where V((b) is the scalar field potential, and the right hand side is supposed to be evaluated at 
the time of entering the long-wavelength regime by a considered mode. (For a recent summary 
of inflation see, for example, [15].) Imagine that the potential V has an inflection point where 
V' = ((f) = 0). Then the predicted amplitude of the generated mode, which happened to enter 
the long-wavelength regime at that moment of time, is infinite. Since for the generation of a 
spectral interval with the Harrison-Zeldovich slope n = 1 one needs the denominator in the r.h.s. 
of (8) to be very close to zero during some interval of time, all the amplitudes in this interval 
of spectrum are predicted by inflationists to be arbitrarily large. Then, according to this logic, 
the relative contribution of gravitational waves becomes negligibly small and formally zero. 

There is absolutely no physical reason for this divergent result in a space-time (gravitational 
pump field) with finite and relatively small curvature, which sets the value of the nominator in 
(8). This result is also in conflict with the finite "temperature" [16] of the De Sitter space-time 
determined by its constant Hubble parameter H. Whatever "particles" are being produced, they 
are supposed to have finite (and small) energy in every frequency interval. However, according 
to (8), even a short interval of the De Sitter evolution is capable of generating an arbitrarily 
large amount of "particles" - density (curvature) perturbations. It was shown [11] that formula 
(8) does not follow from proper evolution and quantum normalization of density perturbations. 
Density perturbations are not simply a scalar test-field, which is normally being considered. 
Density perturbations necessarily involve metric perturbations and do not exist without them. 
At every interval of evolution with the power-law scale factor a(rj), the combined degree of 
freedom satisfies exactly the same equation as eq. (4) for gravitational waves. The associated 
metric (curvature) perturbation h remains constant and small during all the long-wavelength 
regime, similarly to gravitational waves. And every interval of the power-law evolution with 
parameter (3 generates an interval of spectrum with the primordial spectral index n = 2(3 + 5. 
In the correct version of (8), the right hand side must be multiplied by the dimensionless factor 

\J This factor cancels out the zero in the denominator and makes the amplitudes of density 
perturbations finite and of the same order of magnitude as the amplitudes of gravitational waves. 

The real situation on the sky looks even more interesting. It is important to recall [17, 11] 
that a confirmation of the primordial n > 1 ((3 > —2), even at a short spectral interval, will mean 



that the very early Universe was not driven by a scalar field - the cornerstone of inflationary 
considerations - because the n > 1 {(3 > —2) requires the effective equation of state at the initial 
stage to be e + p < 0, but this cannot be accomodated by any scalar field with whichever scalar 
field potential. 

We switch now from cosmology to experimental predictions in terms of laboratory frequencies 
v and intervals of time t (cdt = a(rm)dri). Formula (7) translates into 

/ u \ P+i 

h(u, t) « 1(T 7 cos[2vrz/(t - t v )\ — , (9) 

\VH J 

where vjj = W~ 18 Hz, and t v is a deterministic (not random) function of frequency which does 
not vary significantly on the intervals Au ~ v. We take v = 10 2 Hz as the representative 
frequency for the ground-based laser interferometers. The expected sensitivity of the initial 
instruments at v = 10 2 Hz is h ex = 1CT 21 or better. The theoretical prediction at this frequency, 
following from (9), is h th = 1(T 23 for = -1.8, and h th = 1(T 25 for (3 = -1.9. Therefore, the 
gap between the signal and noise levels is from 2 to 4 orders of magnitude. This gap should be 
covered by a sufficiently long observation time r. The duration r depends on whether the signal 
has any temporal signature known in advance, or not. 

It appears that the periodic structure (9) should survive in the instrumental window of 
sensitivity from u\ (minimal frequency) to v<i (maximal frequency). The mean square value of 
the field in this window is 

/ h 2 {y,t)— = 10- 14 ^-^ / cos 2 [27n^ - t,)]^ +1 du . (10) 

Because of the strong dependence of the integrand on frequency, z^~ 2 6 or ^~ 2 8 , the integral (10) 
is determined by its lower limit. This gives 

ri>2 Ai, / J/-, \ 2/3+2 

/ h 2 (iy,t)— ^10- U — cos 2 [27ri/i(t-*i)] . (11) 

J vi V \Vh) 

The explicit time dependence of the variance of the field, or, in other words, the explicit time 
dependence of the (zero-lag) temporal correlation function of the field, demonstrates that we are 
dealing with a non-stationary process (a consequence of squeezing and severe reduction of the 
phase uncertainty). Apparently, the search through the data should be based on the periodic 
structure at v = v\. 

The response of an instrument to the incoming radiation is s(t) = Fijh 13 where Fij depends 
on the position and orientation of the instrument. The cross correlation of responses from two 
instruments (0\si(t)s2(t)\0) will involve the overlap reduction function [18 - 21], which we assume 
to be not much smaller than 1 [20]. The essential part of the cross correlation will be determined 
by an expression of the same form as (11). 

The signal to noise ratio S/N in the measurement of the amplitude of a signal with no 
specific known features increases as (rz^i) 1 / 4 . If the signal has known features exploited by the 
matched filtering technique, the S/N increases as (r^i) 1 / 2 . The guaranteed law (r^i) 1 / 4 requires 
a reasonably short time r = 10 6 sec to improve the S/N by two orders of magnitude and to reach 
the level of the predicted signal with the extreme spectral index (3 = —1.8. If the law (tz^i) 1 / 2 
can be implemented, the same observation time r = 10 6 sec will allow the registration of the 
signal with the conservative spectral index [3 = —1.9. Even an intermediate law between (rz^i) 1 / 4 
and (rz^i) 1 / 2 may turn out to be sufficient. For the network of ground-based interferometers the 
expected v\ is around 30Hz, but we have used vi = W 2 Hz for a conservative estimate of r. If 
the matched filtering technique can indeed be used, it can prove sufficient to have data from a 
single interferometer. 



For the frequency intervals covered by space intereferometers, solid-state detectors, and 
electromagnetic detectors (see also [22]), the expected results follow from the same formula (9) 
and have been briefly discussed elsewhere [17]. 

In conclusion, the detection of relic (squeezed) gravitational waves may be awaiting only the 
first generation of sensitive instruments and an appropriate data processing strategy. 

Useful discussions with S. Dhurandhar, B. Sathyaprakash, and B. Allen are appreciated. 
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